Susceptibility of a spinon Fermi surface coupled to a U(l) gauge field 
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We study the theory of a U(l) gauge field coupled to a spinon Fermi surface. Recently this model 
has been proposed as a possible description of the organic compound k-(BEDT-TTF)2Cu2(CN)3. 
We calculate the susceptibility of this system and in particular examine the effect of pairing of 
the underlying spin liquid. We show that this proposed theory is consistent with the observed 
susceptibility measurements. 

PACS numbers: 



The organic compounds k-(BEDT-TTF)2X are an in- 
teresting class of materials. Recent experiments have 
shown promise that this compound where the anion, X, is 
Cu2(CN)3 may be the first experimental realized spin liq- 
uid. This material can be described as a nearly isotropic 
effectively two dimensional spin 1/2 triangular lattice at 
half-filling. Experimentally, this material is found to be 
insulating and yet it has no long range magnetic ordering 
down to mK temperatures. Also the static spin suscep- 
tibility remains finite down to the lowest temperatures 
measured. [l| These observations have led to the pro- 
posal that the state may be described by a spinon Fermi 
surface coupled to a [7(1) gauge field. 0, Q The suscep- 
tibility is fit with the high temperature series expansion 
of the spin 1/2 Heisenberg model on a triangular lattice. 
From this fit, the exchange coupling J is found to be 
around 250 K. In addition, the susceptibility is found to 
drop sharply at low temperatures around 10 K before 
saturating to a finite value. [l[ 

Recent measurements of the specific heat have sug- 
gested the existence of a peak in the electronic specific 
heat at around 6 K, once the phonon contribution has 
been subtracted away. [H Led by this discovery, it was 
proposed that the U(l) spin liquid state may have some 
sort of pairing instability. Q Since the specific heat was 
also found to be unaffected by a magnetic field of up to 
8T, conventional singlet pairing is unlikely. The pair- 
ing could, however, be ordinary BCS triplet pairing or a 
new kind of pairing. Recently Lee et al. [B| proposed a 
possible new kind of pairing called "Amperean" pairing. 
Unlike normal BCS pairing across the Fermi surface, this 
pairing is between two spinons on the same side of the 
Fermi surface. In particular, in the Amperean paired 
state, one pairs the spin with momentum Q + p with 
the spin with momentum Q — p where |Q| = fcp and 
|p| small. The Amperean pairing can occur between two 
particles carrying almost parallel momenta due to the 
attractive interaction mediated by the magnetic fluctua- 
tions of the emergent gauge field. As a result the pairs 
carry net momentum 2kp as opposed to in the BCS 
state. In particular, the authors showed that it is possi- 
ble for there to be an instability to this kind of pairing 
for the spinon Fermi surface coupled to a U{1) gauge 
field. They also derived a number of experimental conse- 
quences of this model and show how they could explain 
many of the features seen in the actual experiments on 



k-(BEDT-TTF) 2 Cu 2 (CN)3 Here we calculate the effect 
of pairing on the zero-field spin susceptibility of such a 
system and compare the result to what is experimentally 
seen in this organic compound. 

Starting from a spinon Fermi surface, it is clear that 
at T — the spinons give rise to a Pauli paramagnetic 
term due to the non-zero density of states. Standard 
BCS singlet pairing, however, leads to the reduction of 
this paramagnetism as a gap opens. At first sight, this 
seems to provide a natural explanation of the sharp drop 
in susceptibility below 10 K. However we have already ex- 
cluded BCS singlet pairing because it is inconsistent with 
the observed insensitivity of the specific heat to mag- 
netic field. Both triplet BCS pairing and alternate types 
of pairing such as LOFF and Amperean are consistent 
with the specific heat measurement. However it turns 
out that for such pairing states, the spinon contribution 
to the Pauli paramagnetism is unaffected by the onset of 
pairing, which seems inconsistent with the observed drop 
of susceptibility at low temperatures. In this paper, we 
will show that the drop of susceptibility can be explained 
if the effect of gauge fluctuations is taken into account. 
Before we include the effect of gauge fluctuations, be- 
low we first ignore the gauge fluctuations and explain 
why the onset of pairing does not affect the contribution 
of spinons to the spin susceptibility in the Amperean, 
LOFF and triplet BCS pairing states. 

To see this, we begin with a spinon system with a well 
defined Fermi surface. Applying a magnetic field creates 
two different Fermi seas for the up and down spinons, 
as shown in Fig. [1] First we consider the case of Am- 
perean pairing, where pairing occurs on the same side of 
the Fermi surface. It is possible for both of these spinons 
to lie near the Fermi surface even after the magnetic field 
has been applied (Fig. [1]). This is achieved by pairing 
the spin up spinon with momentum Q + AQ + p with the 
spin down spinon with momentum Q — AQ — p, where 
|p| <C &f an d AQ = (hbH/vf) Q. Moreover, the phase 
space available for p is unchanged with the applied field 
H, as long as the curvature difference between spin up 
and down Fermi surfaces when H ^ can be ignored. 
Thus in this approximation, there is no Zeeman limiting 
field for this pairing. Furthermore, the susceptibility is 
not reduced by pairing because although the opening of 
the pairing gap does smear out the momentum distribu- 
tion rik, it leaves the occupied area of up and down spins 
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FIG. 1: In the Amperean pairing state, the spinons at a and 
b are paired. We see that by choosing AQ + p correctly, both 
a and b lie near their respective Fermi surfaces. 



unchanged. Thus the magnetization remains unchanged 
despite the opening of the pairing gap, and the underly- 
ing spinon paramagnetism is not destroyed. 

We also note that this same argument holds for both 
the LOFF state @, and the BCS triplet state. In the 
LOFF state the pairing is between a spin up on the spin 
up Fermi surface and a spin down on the opposite side 
of the spin down Fermi surface. This state as well only 
smears the momentum distribution leaving the occupied 
area unchanged and again the underlying paramagnetism 
survives. In the particular case of BCS triplet pairing 
where only the equal spin pairings Aj-f and are 
nonzero, it is again clear that the pairing is unaffected 
by magnetic field and thus the magnetization remains 
fixed despite the opening of a pairing gap. In general 
triplet pairing there is also a non-zero pairing between 
the up and down spins. However, due to weak spin-orbit 
coupling the spin quantization axis favors a particular 
direction. When the applied field exceeds some small 
pinning field, the quantization axis rotates in such a way 
that the equal spin pairing description is appropriate and 
the thus the susceptibility remains unchanged despite the 
pairing gap. 

For completeness, we now show how this argument fails 
for the BCS singlet state. After splitting the spin up 
and spin down with a magnetic field, the paired spins 
do not lie on the Fermi surface. If one insists on pairing 
(k, |) with (— k, J.), the gap does not develop at the Fermi 
surface. This is no longer the BCS singlet state, but in 
fact the breeched pairing state, which is not energetically 
favorable for small H. [g] 

Thus we conclude that the naive application of triplet 
or Amperean pairing does not explain the sharp drop in 
the spin susceptibility below 10 K. In order to explain the 
data, we consider the effect of gauge fluctuations in the 
Amperean pairing state. We show that in the problem of 
a spinon Fermi surface coupled to a U(l) gauge field, the 
gauge field fluctuations give rise to a substantial param- 
agnetic contribution to the susceptibility. The onset of 
Amperean pairing driven by the gauge fluctuations in- 
troduces a gap in the gauge field by the Anderson-Higgs 
mechanism and suppress the gauge field contribution to 
the susceptibility. This suppression would account for 
the measured drop in susceptibility that occurs around 
10 K. We note that this mechanism is independent of the 
form of pairing, whether it is triplet or Amperean. 

We begin with the Lagrangian for a 2-D spinon Fermi 



surface system coupled to a U(l) gauge field, 

C = i>*„ (d Q - ia - (x) tp a + —ip* (-iV - a) (1) 

2m 

We have dropped the gauge field kinetic energy term be- 
cause its strength is inversely proportional to the charge 
gap. ip a is the spinon field and the gauge field is 
a = (ao, a), fx is the chemical potential. We work in 
Coulomb gauge V • a = 0. 

We begin by calculating the random-phase approxi- 
mation (RPA) for this model. The use of the RPA can 
be justified by the standard 1/N expansion. 0, [l(| El 
We need to calculate the bare spinon polarization bubble 
shown in Fig. [4] which generates the gauge propagator. 
Working with this bare spinon bubble, i.e. not dressed 
by further fields inside the spinon loop, is equivalent to 
working to lowest order in the 1/N expansion. The scalar 
i.e. longitudinal part of the gauge propagator is related 
to the density-density response and does not show sin- 
gular behavior for small q and to. In other words the 
scalar part is screened out by spinon density fluctuations 
and we can focus on just the transverse part of the gauge 
field. Because the gauge field is now purely transverse, 
the spinon-gauge field vertex carries a vector index. 

The gauge propagator is generated by a sum of spinon 
loops carrying a given spin. We define the bare spinon 
polarization bubble to be 
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fco 



k + | ) G a ( I, 



k x q 



( 3 ) 

with q = (qo, q) and k similary defined. We show the 
details of this calculation in the appendix, but in the end 
we find that 
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ky and Vp are the Fermi 



where Yd = tt; — and 7 = 
momentum and velocity respectively for a spinon with 
spin a in an applied magnetic field. We will see that we 
need to keep the curvature of the Fermi surface so we use 
that kp — y/kp ± 2mfiBH, where a = (+, — ) for the up 
and down spins respectively. 

In order to calculate the susceptibility, we calculate the 
free energy of the Lagrangian from Eq. [T] We first note 
that because of the vector nature of the gauge propaga- 
tor, there are no tadpole diagrams. The diagram we con- 
sider then is the standard RPA which is a closed string of 
bubbles. This is equivalent to calculating the free energy 
by integrating out the spinon fields and obtaining an ef- 
fective action for the gauge field. The effective action is 
5(a) = }~2 q n(g)aja q with 11(g) from Eq. H Thus the 
partition function is 



Z = e-? F = / Dae'*) 



(5) 
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where F is the thermodynamic potential. Performing the 
functional integral, we find 



F 



-T 

R ^ 



lnll(g). 
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We need to calculate the susceptibility \ — — fjp-|//=o- 
The spinon bubble in Eq. [2] only depends on H through 
the parameters Vp and fcp of W a (q) . Inserting Eq. [H into 
Eq. [21 we then Taylor expand to get 



n(<z) = n (q) + A(q)H 2 + 0(H 3 



where 
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and 
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Plugging this expansion into the definition of the sus- 
ceptibility, we find that the correction to the zero field 
susceptibility due to this RPA diagram is 
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From the above functional forms of A(q) and Ilo(g), it 
is clear that Ax > 0. Thus the RPA correction to the 
gauge field gives rise to an additional paramagnetic con- 
tribution to the total susceptibility. When the gauge field 
is gapped after pairing, this correction is reduced. This 
is consistent with the measured drop that is seen exper- 
imentally. 

An alternate way of calculating the spin susceptibil- 
ity would be to dress the spin-spin correlation function 
with the gauge field propagator from Eq. 2J From this 
one could calculate the susceptibility explicitly. These 
two methods are the same because given the free energy 
RPA diagrams, taking each derivative with respect to 
magnetic field is equivalent to adding a spin-flip vertex. 
By writing down all the topologically inequivalent ways 
to add two spin-flip vertices, one generates the diagrams 
of the spin-spin correlation function. Kim et. al. cal- 
culated the density-density correlation function for small 
q and found that the singular portions coming the self- 
energy and vertex corrections cancel and that there are 
only analytic corrections. [12j This result applies equally 
to the spin-spin correlation function. Our calculation cor- 
responds to calculating the numerical value of the non- 
singular part of this correction. 

We now proceed to calculate the numerical value of the 
shift in the susceptibility in this RPA approximation for 
the gauge propagator. From the derived forms of Ha(q) 
and A(q), we see that the integrand of Eq. [TO]increases as 
q increases. Thus to calculate a numerical value we need 
to introduce a cutoff. In the derivation of the polarization 
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FIG. 2: Ax in units of xo as a function of the pairing gap 
Eg in units of ep. The solid, dotted, and dashed lines are for 
cutoffs 1,1.5 and 2&f. 



bubble in Eq. [4] we assumed that q < fcp. Thus we let 
the energy integration go from zero to infinity and cutoff 
the momentum integral at something of order kp . In the 
standard calculation of the density-density polarization 
bubble for fermions in two-dimensions 13], one finds that 
the polarization bubbles dies off sharply as a square root 
for q > 2fcp and is relatively flat inside the Fermi surface. 
Thus in order to compare to experiment, we take the 
cutoff to be 2kp. We write the paramagnetic correction 
at the RPA level in terms of the 2-D Pauli paramagnetic 
susceptibility for a spin 1/2 free fermion system, xo = 
m/V, where m is the effective mass of the spinons. For 
this cutoff, the extra paramagnetic contribution is A% ~ 
0.42 xo- This number, however, is strongly dependent on 
the cutoff as seen in Fig. 

We now introduce the gapping of the gauge field due 
to spinon pairing. Calling the gap energy E g , once the 
spinons begin to pair, 11(g) is shifted by this energy and 
as a result the denominator of the integrand of Eq. [TU] 
becomes IIo (q) + E g . As the temperature decreases E g 
rises and the paramagnetic correction falls. In Fig. [21 we 
plot Ax (-E s ), the extra paramagnetic contribution due to 
these RPA diagrams, for different values of the cutoff. 

We now consider the effect of the short range interac- 
tion on this calculation. We work perturbatively in U, 
the strength of this interaction. We need to consider the 
free energy diagrams that contain both the gauge field 
and the short range interaction. We continue to work us- 
ing the RPA. To lowest order in U there are four diagrams 
that contribute to the propagator by adding corrections 
to the spinon bubble 11(g). These diagrams are shown in 
Fig. [3] The top two diagrams are the standard Hartree 
and Fock corrections to the spinon propagator. These 
two are clearly related by an exchange. One can also see 
that the bottom two diagrams in Fig. [3] are also related 
to each other via an exchange. 

First we consider the bottom two diagrams. We denote 
the correction to the total polarization for the vertex cor- 
rection diagram (Fig. [3b) as II C and for its associated ex- 
change diagram (Fig. [3JI) as lid- In order to calculate the 
free energy, we close these diagrams with the propagator 
calculated from the bare spinon bubble D(q) = 11(g) -1 . 
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FIG. 3: Lowest order corrections in the short range interac- 
tion, U, to the gauge propagator. 



We find that 

U c ( q ) = J dk dk' G{k)G{k + q)G{k')G{k' + q) x 
[k-k'-(k-q)(k'.q)l 



'-U(k-k'), (11) 



dk dk' G(k)G(k + q)G(k')G(k' + q) x 



[k-k'-(k.q)(k'.q)] 



U(q). (12) 



Assuming, for simplicity, that U(q) is independent of q, 
it is clear that for both diagrams the integral over #kq, 
the angle between k and q, is odd and thus gives zero. 
This is because unlike in Eq. [4] there is now only one 
power of k from the vertex. We note that even without 
the assumption of U being momentum independent lid 
is zero, but that the assumption is necessary for Il c . 

The short range interaction can also dress the spinon 
propagator with Hartree-Fock corrections. Instead of 
evaluating these diagrams explicitly, we use the Dyson 
equation to calculate the dressed spinon propagator for 
a given spin, G a . The Hartree term, diagram Fig. [3Jd, 
gives a self energy Un where n is the total density of 
spinons. The Fock term, diagram Fig. [3^, gives a self 
energy U n a . Thus we have 



G„ = 



1 



(13) 



The total spinon density does not depend on magnetic 
fied, thus the Hartree term contributes a field indepen- 
dent shift to the chemical potential which does not effect 
the susceptibility and thus is dropped. Applying a mag- 
netic field does shift the density of a particular spin to 
n a {H) = n CT (0) + j^Sfi a: where Sfi a = ±.2^H for up 

and down spins respectively and ^ = ^ is the standard 
2-D density of states. Thus, 



G„ — 



1 



(e k - (jjl + J/v) - U(n a (0) + |J<W)) 



Again, dropping a shift in the chemical potential, this 
dressed spinon propagator is the same as the undressed 
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FIG. 4: The one loop corrections to the bare gauge propaga- 
tor. 



one except that when a magnetic field H is applied it 
responds to an effective field H = H (l + ^jf)- 

The gauge-field correction to the susceptibility includ- 
ing the short-range interaction to lowest order x is thus 




(15) 
(16) 



where \ is the paramagnetic correction from Eq. 1101 
Thus for U > 0, we see that including the short range 
interaction diagrams enhances the extra paramagnetic 
term and thus leads to a larger drop in the susceptibility 
as the pairing gap opens. 

In conclusion, we have shown that the proposed 
Amperean pairing of a spinon Fermi surface coupled 
to a U(l) gauge field is consistent with the experi- 
ments performed on the candidate material k-(BEDT- 
TTF) 2 Cu2(CN) 3 . In particular, the unconventional pair- 
ing of spinons on the same side of the Fermi surface allows 
for a non-zero T = susceptibility despite the opening 
of a gap. Also by calculating the effect of the gauge 
field on the paramagnetic susceptibility, we found a drop 
in the susceptibility as the gauge field becomes gapped 
due to pairing. This is consistent with the drop seen in 
the experiments. However, since the contribution to the 
susceptibility comes from the gauge field carrying large 
momentum, the result is sensitive to the cut-off and the 
exact numerical factor cannot be trusted. Our goal is 
rather to show as a matter of principle that there is a 
large paramagnetic contribution to the susceptibility that 
is suppressed by the onset of pairing. 

We thank Y. Nakazawa and K. Kanoda for sharing 
their data prior to publication. P.A.L acknowledges sup- 
port by NSF grant DMR-0517222. 



APPENDIX A: DERIVATION OF 11(g) 

In this appendix, we present the details of the calcula- 
tion of the gauge field propagator generated by the bare 
spinon bubble, Eq. [U We define U(q) = IL a (q) + H b (q), 
where the two one-loop diagrams that contribute are 
shown in Fig. [4] 
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We start with the diagram in Fig. [4^. Thus 



kg 
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where the extra factor of two came from the spin sum- 
mation. Just as in the standard calculation of the polar- 
ization bubble we do the Matsubara sum first and get 



n o (g) = -2 
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where = efe — /i, efe being the spinon dispersion relation 
and /x the Fermi energy. If |q| < k-p then the difference 
of the two Fermi functions is only nonzero for a region 
of length qcos{9) around kp. Assuming that q is much 
smaller than k-p, the k integration is just replaced by 
k = k-p. We now calculate the real and imaginary parts 
of Tl(q) and find 
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-k 3 
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2it 2 m 2 J q 2 + vlq 2 cos 2 (6») 



Doing these integrals, we find that the imaginary part is 
zero and that the real part gives, 
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where A± = ±5^ + fm' We Taylor expand each distri- 
bution function in A±. We work to third order in A± in 
order to generate all the terms up to order q 2 . Dropping 
all higher powers of q, we find that in the static limit, 
90 = 0, 
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Which gives 



n„(<7) = — 
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We now compute the effect of the diagram shown in 
Fig. IDd. The spinon loop here just gives a density. So 
that this diagram evaluates to 



n 6 (g) = — (-n) = 



(A9) 



Putting Eq. 
that 



5J Eq. IA8I and Eq. IA9I together, we get 



n(g) 
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v/i; 2 q 2 + ql + q 
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This gives the dominate correction for small q for the 
non-static polarization bubble. For the static part of the 
polarization bubble, qo = 0, this term vanishes, so we 
have to relax the approximation to the difference of the 
Fermi functions in order to get the q dependence in the 
static case. We now write 



4±h = £ k + A d 
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where 7 



and Xd 



12m7r * 



Taking the limit go «C vpq, we recover the result from 
earlier papers 11(g) = + Xdg 2 - 0, E3] Note that the 
actually value of Xd differs from [ijj by a factor of two 
because in that paper, there is a factor of two error in the 
form of the 2-D Landau diamagnetic susceptibility that 
gives rise to Xd- 
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